Multi-band Gutzwiller wave functions for itinerant ferromagnetis m0 
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Multi-band Gutzwiller-correlated wave functions reconcile the contrasting concepts of itinerant 
band electrons versus electrons localized in partially filled atomic shells. The approximate evaluation 
of these variational ground states becomes exact in the limit of large coordination number. The 
result allows the identification of quasi-particle band structures for correlated electron systems. As 
a first application, we summarize a study of itinerant ferromagnetism in a two-band model, thereby 
elucidating the co-operation of the Coulomb repulsion and the Hund's-rule exchange. Then, we 
present results of calculations for ferromagnetic nickel, using a realistic 18 spin-orbital basis of 4s, 
4p and 3d valence electrons. Good agreement with the experimental ground-state properties of 
nickel is obtained. In particular, the quasi-particle energy bands agree much better with the photo- 
emission and Fermi surface data than the band structure obtained from spin-density functional 
theory. Finally, we present results for the variational spinwave dispersion for our two-band model. 



I. INTRODUCTION 

How well do we understand the itinerant ferromag- 
netism of iron, nickel, and other transition metals and 
their compounds? 

More than 50 years ago two basically different scenarios 
had emerged from early quantum-mechanical considera- 
tions on electrons in metals with partly filled d bands. In 
scenario I, suggested by Slater jij and Stoner |Q, band 
theory alone was proposed to account for itinerant fer- 
romagnetism: due to the Pauli principle, electrons with 
parallel spins cannot come arbitrarily close to each other 
("Pauli" or "exchange hole"), and, thus, a ferromag- 
netic alignment of the electron spins reduces the total 
Coulomb energy with respect to the paramagnetic situ- 
ation ("exchange energy"). This scenario, equivalent to 
Hartree-Fock theory, was criticized by van Vleck ||. In 
an independent-electron theory electrons are distributed 
statistically over the lattice which implies strong charge 
fluctuations on the transition metal atoms. 

Consequently, as a scenario II, van Vleck || empha- 
sized the importance of electron correlations: due to the 
strong electron-electron interactions, charge fluctuations 
in the atomic d shells are strongly suppressed ("mini- 
mum polarity model"), and, thus, atomic rather than 
band aspects are decisive for itinerant ferromagnetism. 
In scenario II, the atomic magnetic moments arise due to 
the local Coulomb interaction (in particular, Hund's-rule 
couplings), and they may align because of the electrons' 
motion through the crystal. 

In principle, such a dispute can be resolved in nat- 



ural sciences. The corresponding theories have to be 
worked out in detail, and their results and predictions 
have to be compared to experiments. This was indeed 
done for scenario I [jjj. The (spin-)density functional the- 
ory is a refined band theory which describes simple met- 
als with considerable success. Unfortunately, progress 
for scenario II was much slower. It calls for a theory of 
correlated electrons, i.e., a genuine many-body problem 
has to be solved. It was only recently that reliable the- 
oretical tools became available which allow to elucidate 
scenario II in more detail §,|,0J|§. 

A first step in this direction was the formulation of 
appropriate model Hamiltonians which allowed to dis- 
cuss matters concisely. In 1963/1964 Gutzwiller |Io| , pT| , 
Hubbard {T^f , and Kanamori []13| independently intro- 
duced the Hubbard model which since then has become 
the standard model for correlated electrons on a lattice. 
This model covers both aspects of d electrons on a lattice: 
they can move through the crystal, and they strongly in- 
teract when they sit on the same lattice site. The model 
is discussed in more detail in Sec. ||. 

Even nowadays, it is impossible to calculate exact 
ground-state properties of such a model in three dimen- 
sions. In 1963/1964 Gutzwiller introduced a trial state to 
examine variationally the possibility of ferromagnetism in 
such a model lOpTj]. His wave function covers both lim- 



its of weak and strong correlations and should, therefore, 
be suitable to provide qualitative insights into the mag- 
netic phase diagram of the Hubbard model. Gutzwiller- 
correlated wave functions for multi-band Hubbard mod- 



els are defined and analyzed in Sec. Ill 



*Dedicated to Martin C. Gutzwiller on the occasion of his 75th birthday. 
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As envisaged by van Vleck |J, Gutzwiller [|0| found 
that ferromagnetism is non-generic in the one-band Hub- 
bard model, in contrast to the predictions of the Stoner 
theory. The reason is simple: the Gutzwiller wave func- 
tion not only describes the exchange hole between elec- 
trons of the same spin species but also the "correlation 
hole" between electrons of different spin species which 
forms when the electron-electron interaction is taken into 
account consistently. Therefore, the energy difference 
between the correlated paramagnetic and ferromagnetic 
states is very small, and paramagnetism wins in almost 
all one-band cases. In the multi-band situation |ll|] , itin- 
erant ferromagnetism is favored by two concomitant ef- 
fects: (i) the suppression of atomic charge fluctuations 
such that, (ii), the local exchange interactions give rise 
to a local magnetic moment as in isolated atoms (Hund's 
rule). Hence, as proposed by van Vleck, small charge fluc- 
tuations and local exchange interactions generate large 
local moments which are present in both the paramag- 
netic and ferromagnetic phases. At low temperatures the 
electrons' motion through the crystal may then lead to 
long-range order of these pre- formed local moments. In 
this way, Gutzwiller's approach unified and substantiated 
van Vleck's early ideas on a correlated-electron theory of 
itinerant ferromagnetism. 

Unfortunately, the evaluation of multi-band Gutzwiller 
wave functions itself poses a most difficult many-body 
problem. Perturbative treatments p||l5| ] are constrained 
to small to moderate interaction strengths. The region 
of strong correlations could only be addressed within the 
"Gutzwiller approximation" j^|ll],[l6| and its various ex- 
tensions p7j , ^8| . Thus, an application to ferromagnetism 
of real materials was not in sight, and Gutzwiller's qual- 
itative findings were not appreciated much in the sixties. 

About 25 years later, a major technical advancement 
opened a new way to analyze Gutzwiller-correlated wave 
functions. The one-band Gutzwiller wave function was 
evaluated exactly in one dimension |]l9|j20|| , and the 
Gutzwiller approximation was found to become exact for 
the one-band Gutzwiller wave function in the limit of in- 
finite spatial dimensions, d — ► oo p^B^B^. Based on 
a diagrammatical approach, Gebhard |23| developed a 
compact formalism which allows to calculate the varia- 
tional ground-state energy in infinite dimensions without 
the evaluation of a single diagram. He found that \/d cor- 
rections are quantitatively small for the Gutzwiller wave 
function, i.e., the result in infinite dimensions is a reliable 
approximation to three dimensions. 

Recently, Gebhard's approach was generalized by us 
to the case of multi-band Gutzwiller wave functions J9|. 
Thereby, earlier results by Biinemann and Weber [B4|, 
based on a generic extension of the Gutzwiller approxi- 
mation Bah were found to become exact in infinite di- 
mensions ]26|. Half a century after van Vleck's call for 
new many-particle methods and more than thirty years 
after Gutzwiller's first paper on the subject, techniques 



are at hand to carry out the program outlined above as 
scenario II. 



Our results on a two-band toy model in Sect. IV con- 
firm Gutzwiller's insights on the fundamental require- 
ments for itinerant ferromagnetism. Substantial atomic 
Coulomb interactions suppress local charge fluctuations 
and thereby allow the generation of local moments due to 
the atomic Hund's-rule exchange couplings. Our quan- 
titative results from a full-scale calculation for nickel in 
Sect. |y| confirm these qualitative findings. The quasi- 
particle properties of this material (e.g., exchange split- 
tings, density of states, topology of the Fermi surface) are 
reproduced quantitatively from our Gutzwiller-correlated 
wave functions. In particular, strong electron correla- 
tions are required to describe nickel appropriately. 



In Sect. VI, we briefly discuss a new method to cal- 



culate spinwave spectra of correlated electron systems 
based on the Gutzwiller variational scheme ^7j. In 
this way, the dispersion relation of the fundamental low- 
energy excitations can be derived consistently. Albeit 
the description is based on itinerant electrons, the re- 
sults for strong ferromagnetism in a generic two-band 
model resemble those of a Heisenberg model for localized 
spins. Again, the Gutzwiller approach provides a com- 
prehensive scheme for a unified description of localized 
and itinerant aspects of transition metals. 

In this brief review we skip all technical details and 
refer the reader to the original literature for further in- 
formation. 



II. HAMILTON OPERATOR 

Our multi-band Hubbard model JT3] is defined by the 
Hamiltonian 



H 



2^ tif^c 



Hi [a t = Hi + H a 



(1) 



Here, c^ a creates an electron with combined spin-orbit 
index a = 1,...,2N (N = 5 for 3d electrons) at the 
lattice site i of a solid. 

The most general case is treated in Ref. ||. In this 
work we assume for simplicity that different types of 
orbitals belong to different representations of the point 
group of the respective atomic state (e.g., s, p, d(e g ), 
d(t2g))- In this case, different types of orbitals do not 
mix locally, and, thus, the crystal field is of the from 



t 



c u er,<T 



Consequently, we use normalized single- 
particle product states |$q) which respect the symmetry 
of the lattice, i.e., 



(2) 



We further assume that the local interaction is site- 
independent 
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(3) 



This term represents all possible local Coulomb interac- 
tions. 

As our basis for the atomic problem we choose the con- 
figuration states 



10 = |0"1,0"2, 



c *>i c *>3 ' ' ' Kaccum) (m < <r 2 < 



(4) 



which are the "Slater determinants" in atomic physics. 
The diagonalization of the Hamiltonian i^ ;a t is a stan- 
dard exercise ^8|. The eigenstates |F) obey 



\T) =^^1) , 



(5) 



where Tj -p are the elements of the unitary matrix which 
diagonalizes the atomic Hamiltonian matrix with entries 
(I\H iM \r). Then, 



r 

m r = |r)(T|. 



(6) 
(7) 



The atomic properties, i.e., eigenenergies Ep, eigen- 
states |r), and matrix elements Tjr, are essential in- 
gredients of our solid-state theory. 



III. MULTI-BAND GUTZWILLER WAVE 
FUNCTIONS 

Gutzwillcr-corrclatcd wave functions are written as 
a many-particle correlator Pq acting on a normalized 
single-particle product state |$o)> 



(8) 



The single-particle wave function |$o) contains many 
configurations which are energetically unfavorable with 
respect to the atomic interactions. Hence, the correla- 
tor Pq is chosen to suppress the weight of these con- 
figurations to minimize the total energy in ([!]). In the 
limit of strong correlations the Gutzwiller correlator Pq 
should project onto atomic eigenstates. Therefore, the 
proper multi-band Gutzwiller wave function with atomic 
correlations reads 

i 

ho = II \p T = I] I 1 + ( A -r - 1) fci;r] (9) 
r r 

= 1 + ^ (A i;r - 1) m i;r . 
r 

The 2 2N variational parameters Ai ; r per site are real, 
positive numbers. For Ai ; r and all other A; ; r = 



all atomic configurations at site i but \Yq) are removed 
from |$o)- Therefore, by construction, I^g) covers both 
limits of weak and strong coupling. In this way it in- 
corporates both itinerant and local aspects of correlated 
electrons in narrow-band systems. 

The class of Gutzwiller-correlated wave functions as 
specified in (j^) was evaluated exactly in the limit of infi- 
nite dimensions in Ref. Q . The expectation value of the 
Hamiltonian reads 



(H) 



(*g\H\* g ) 
(*g|*g) 

E / 



i\a 



f- ^2 E r m v j 
i;r 



(10) 



Here, = ($o|^i;o-|^o) is the local particle density in 
|$o)- The local q- factors are given by ||] 



n°(l-n°) ^ 

(TV (7 / p p/ 



mrmr / 



E 



fa fa 



V TO (/'ua)< T r,(7u CT )Vu<.) ) rT+ >J ,T I>r( , (11) 

where mj.j (mj r ) is the probability to find the config- 
uration | J) (the atomic eigenstate on site i in the 
single-particle product state l^o)- The fermionic sign 
function 



/i = (7U CT |c+|J) 



(12) 



gives a minus (plus) sign if it takes an odd (even) num- 
ber of anticommutations to shift the operator to its 
proper place in the sequence of electron creation opera- 
tors in \I U a). 

Eqs. (0) and ( fill ) show that we may replace the orig- 
inal variational parameters A; ; r by their physical coun- 
terparts, the atomic occupancies m !; r. They are related 
by the simple equation bl 



rrii-r = Xl r ml 



i-,r ■ 



(13) 



The probability for an empty site (|/| — 0) is obtained 
from the completeness condition, 



E m ^ r 

r(|r|>i) 



(14) 



The probabilities for a singly occupied site (|/| = 1) are 
given by 



rrii-o 



n° 



I(m>2)(ae/) 



K 



EE 



(15a) 
(15b) 
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The parameters m i; g and m^ a must not be varied in- 
dependently. All quantities in ( |To| ) are now expressed 
in terms of the atomic multi-particle occupancies rrii-r 
(|r| > 2), the local densities n®. a , and further variational 
parameters in |$o)- 

It is seen that the variational ground-state energy can 
be cast into the form of the expectation value of an effec- 
tive single-particle Hamiltonian with rcnormalizcd elec- 
tron transfer amplitudes t\ 



l C7 ,CT 



(16) 

Therefore, |$o) is the ground state of H c g whose pa- 
rameters have to be determined self-consistently from 
the minimization of ($o|if e g |$o) with respect to m^r 
and n®. a . For the optimum set of parameters, de- 
fines a band structure for correlated electrons. Similar 
to density-functional theory, this interpretation of our 
ground-state results opens the way to detailed compar- 
isons with experimental results. 

Applications are discussed in the next two sections. 



IV. RESULTS FOR A GENERIC TWO-BAND 
MODEL 

The atomic Hamiltonian for a two-band model (b — 
1, 2) can be cast into the form 

Hi-^t = U 2J ?i>b,inb,i + U' 2J hi^h 2 ,a> - J } j hi^h 2 ,<, 

b cr.fr' fJ 



2,-cr l,-cr u 2,cr 



(17) 



+ Jc {^1,1^1,1^2,1^2,1 ^2, 1^2,1^1,1^1,1 



For two orbitals, H at exhausts all possible two-body in- 
teraction terms. 

We assume that the model describes two degenerate 
d(e g ) orbitals which leads to the following restrictions 
enforced by the cubic symmetry p8[ : (i) J = Jq, and 
(ii) U — U' = 2 J. Therefore, there are two independent 
Coulomb parameters, the local Coulomb repulsion U and 
the local exchange coupling J. Since U,U' — O(10 eV) 
are of the same order of magnitude the relation (ii) shows 
that J = 0(1 eV) is of the typical strength of the atomic 
Hund's rule couplings. 

For the one-particle part Hi we use an orthogonal 
tight-binding Hamiltonian with first and second nearest 
neighbor hopping matrix elements. We apply the two- 
center approximation for the hopping matrix elements 
and exclude any spin-flip hopping. The resulting band 



structure (width W = 6.6 eV) is discussed in detail in 
Ref. @. 

In the following we concentrate on two band- 
fillings < n < 4, the number of electrons per lat- 
tice site. Alternatively, we use < n„ = n/4 < 1, 
the particle density per band and per spin direction. 
For n a = 0.29, the non-interacting density of states has 
a maximum which is most favorable for ferromagnetism 
according to the Stoner (i.e, Hartree-Fock) theory. More- 
over, we study n a = 0.35, where the density-of-states has 
a positive curvature. 
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FIG. 1. Phase diagram as a function of U and J for the 
Hartree-Fock-Stoner theory (HF) and the Gutzwiller wave 
function (GW) for (a) n/4 = 0.29 and (b) n/4 = 0.35; PM: 
paramagnet, FM: ferromagnet. 

In Fig. [j] we display the J-U phase diagram for 
both fillings. It shows that Hartree-Fock theory always 
predicts a ferromagnetic instability. In contrast, the 
correlated-electron approach strongly supports the ideas 
of van Vleck || and Gutzwiller ]ll|j : (i) a substantial 
on-site exchange J is required for the occurrence of fer- 
romagnetism if, (ii), realistic Coulomb repulsions U are 
assumed. In 1964 Gutzwiller wrote ". . . one may hope to 
show some day that the terms (ii) can never induce fer- 
romagnetism at all by themselves . . . " . Fig. [j] strongly 
supports Gutzwiller 's ideas on the importance of the lo- 
cal Hund's-rule couplings. At the same time the compar- 
ison of Figs, [l]a and [j]b shows the importance of band- 
structure effects which are the basis of the Stoner the- 
ory. The ferromagnetic phase in the U-J phase diagram 
is much bigger when the density of states at the Fermi 
energy is large. Therefore, the Stoner mechanism for fer- 
romagnetism is well taken into account by Gutzwiller's 
correlated-electron approach. 
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FIG. 2. Condensation energy as a function of U for 
J = 0.2(7 for the Hartree-Fock theory (HF) and the 
Gutzwiller wave function (GW) for n/4 = 0.29 (full lines) 
and n/4 = 0.35 (dashed lines). 

In Fig. |^, we display the energy differences between 
the paramagnetic and ferromagnetic ground states ( "con- 
densation energy", -E C ond) as a function of the interac- 
tion strength for J = Q.2U. This quantity should be 
of the order of the Curie temperature which is in the 
range of 100 K — 1000 K in real materials. The Hartree- 
Fock-Stoner theory yields such small condensation ener- 
gies only in the range of U » 4eV; for larger U, E con d 
is of order U . In any case, the interaction parameter U 
has to be tuned very precisely to give condensation en- 
ergies which concur with experimental Curie tempera- 
tures [Q. In contrast, for the Gutzwiller-correlated wave 
function, we find relatively small condensation energies 
-Econd — 0.5 • 10 3 K even for interaction values as large 
as twice the bandwidth (U m 12 eV). Moreover, the de- 
pendence of the condensation energy on U is rather weak 
such that uncertainties in U do not drastically influence 
the estimates for the Curie temperature. 

These toy-model studies show that the Gutzwiller vari- 
ational approach contains all the fundamental elements 
for a successful theory of itinerant ferromagnetism. In 
the next section we show how to quantify this statement 
for ferromagnetic nickel. 

V. CORRELATED BAND-STRUCTURE OF 
NICKEL 

A. Discrepancies between experiment and SDFT 



the Fermi energy (3.0 versus 1.9 states/(eV atom)), which 
indicates a quasi-particle mass enhancement by a factor 
of approximately 1.6. Here, the Sommcrfcld formula is 
used to convert the specific heat data; the theoretical 
value follows directly from the quasi-particle band struc- 
ture. Furthermore, very detailed photo-emission studies 
at symmetry points and along symmetry lines of the Bril- 
louin zone show discrepancies to SDFT results for indi- 
vidual band-state energies which are of similar magnitude 
as seen in the overall d bandwidth. 

The studies also revealed even bigger discrepancies 
in the exchange splittings of majority spin and mi- 
nority spin bands. The SDFT results give a rather 
isotropic exchange splitting of about 600 meV . 
In contrast, the photo-emission data show small and 
highly anisotropic exchange splittings between 160 meV 
for pure d(e g ) states such as X 2 and 330 meV for pure 
d{t2 g ) states, the latter value estimated from the ex- 
change splitting of A3 states along V to L (|^J3^]. The 
much larger and much too isotropic exchange splitting of 
the SDFT results has further consequences. 

1. The experimental magnetic moment of the strong 
ferromagnet Ni is /i = 0.61/ib; yet of relevance is its 
spin-only part /i S pm-oniy = 0.55ub (34|]. The SDFT 
result is /i S pin-oniy = 0.59/^b ||, an overestimate 
related to the too large exchange splitting. 

2. the X2 state of the minority spin bands lies below 
E-p |3q| , whereas all SDFT results predict it to lie 
above the Fermi level ||,[36],[37]] . As a consequence, 
the SDFT Fermi surface exhibits two hole ellipsoids 
around the X point of the Brillouin zone while in 
the de-Haas-van-Alphen experiments only one el- 
lipsoid has been found p3|]3"s| ] . 

3. The strong t2 g -e g anisotropy is also reflected in the 
total d hole spin density, i.e., in the observation 
that the d-hole part of the Ni magnetic moment 
has 81% d(t2 9 ) and 19% d(e g ) character Q, while 
the SDFT results give a ratio of 74% to 26% |§. 

In the late 70's and early 80's various authors have 
investigated in how far many-body effects improve the 
agreement between theory and experiment, see, e.g., 

Refs. mm. For example, Cooke et al. [Q introduced 

an anisotropic exchange splitting as a fit parameter. 



Of all the iron group magnetic metals, nickel is the 
most celebrated case of discrepancies between the results 
from experiment and from SDFT. From very early on, 
the photo-emission data have indicated that the width 
of the occupied part of the d bands is approximately 
W* cc = 3.3 eV || whereas all SDFT results yield val- 
ues of W * cc SDFT = 4.5 eV or large r P , p9| . Similarly, the 
low temperature specific heat data |3C| ] give a much larger 
value of N*(Ep), the quasi-particle density of states at 



B. Parameterization and minimization 

For the single-particle Hamiltonian Hi in (Q) we use 
a nine-orbital s, p, d basis. The hopping matrix ele- 
ments t1'° are determined from a least-square fit to en- 
ergy bands obtained from a density-functional-theory cal- 
culation for non-magnetic nickel. The root-mean-square 
deviation of the d-band energies is about 5 meV. Details 
of this calculation will be published elsewhere E2| . 
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The atomic Hamiltonian H at in (||) is restricted to 
only d orbitals. As a consequence, there are three in- 
dependent interaction parameters (we use the spherical- 
atom approximation), which can be expressed either as 
Slater integrals F k (fc = 0,2,4) or as Racah parame- 
ters A, B, C j28|. We neglect any cubic crystal-field cor- 
rections to the interaction parameters and, furthermore, 
put C/B = 4.5, as suggested by the results of ligand- 
ficld theory @ . We find that A » 10 eV and C « 0.1 eV 
are the choice which best reproduces the experimental 
results; see below. 

The number of multi-electron states |T) is 2 10 — 11; 
because of the cubic site symmetry, the number of in- 
dependent (internal) variational parameters mr reduces 
to approximately 200 for paramagnetic and to approx- 
imately 400 for ferromagnetic solutions. For the mag- 
netic solutions, there exist a certain number of (exter- 
nal) variational parameters A a in the one-particle wave 
function |$o)- in particular, the A Q influence the spin- 
orbital densities n\. b a (b — d,p,s). Among the A Q are 
the "exchange splittings" of the individual orbitals (b, a) 
and (b, — a), i.e., these parameters determine the degree 
of magnetization in |3>o). Most important are the ex- 
change splittings for d(e g ) and dfag) orbitals, of minor 
importance are those for s and p orbitals. Furthermore, 
the d(e g ) vs. d(t2g) crystal-field splitting has to be in- 
cluded in the set of A a . Technically, each |$o) , which is a 
functional of the external variational parameters A Q , de- 
fines a band structure and, consequently, a Fermi surface. 
Therefore, calculating the particle densities n°. fe a and 

the expectation values ($o|-ffcff |^o) implies momentum- 
space integrations up to the respective Fermi surface. It 
should be noted that the introduction of the variational 
parameters A a causes a charge flow between the d, and 
s and p states which has to be compensated by using 
appropriate "chemical potentials" . 

The optimum |$ ) for a given set of {A Q }, |$Q Pt ) a is 
the ground state of the effective single-particle Hamil- 
tonian (|l^) with renormalized hopping matrix elements 
tjj which, through the q factors, are functions of the in- 
ternal variational parameters mr- Thus, |$Q Pt ) Q can be 
obtained from a self-consistent procedure, starting with 
q a = 1, i.e., with |<I>q e ) a . After minimization of the 
internal variational parameters, the resulting values for 
q a define a new effective Hamiltonian from which a new 
|$o)q is constructed. In this step, the Fermi surface inte- 
grations have to be repeated. Self-consistency is usually 
reached rather quickly, i.e., |$o Pt ) Q is found after three 
to five iterations. 

The global minimum, I^q 15 ') global is found by a search 
through the A^ Q -dimensional space of the external vari- 
ational parameters A a . This search can be sped up by 
first optimizing with respect to the most important ex- 
ternal variational parameter which is the isotropic ex- 
change splitting Ad- In a second step, the anisotropy of 



the exchange splitting is investigated, i.e., we introduce 
A eg and A t2g , keeping the average A<j close to the value 
of A° pt . The searches for A° pt and for A°p' and A? pt 
can be carried out starting with |$Q arc ). Only then the 
self-consistency procedure for |3>Q Pt ) has to be launched. 
Typical energy gains are (in meV): 

^barc _ £;baro( A opt) _ jq _ ^ ^g^ 
£ ,ba rc(A opt ) _ £ bare (A opt 5 A opt } (lgb) 
E b-e (A opt ) A opt } _ £opt (A op t) A opt) , 5 _ 10 . (l 8c) 

The energy gains from the variations of A s and A p are 
of the order of 0.1 meV. 

At present, our description does not contain the spin- 
orbit coupling which could be relevant in some transition 
metals. This interaction contributes to the orbital mo- 
ment ^(orbit = 0.06/ie to the total magnetic moment. The 
inclusion of the spin-orbit coupling poses no principle dif- 
ficulties ||] but it considerably enhances the numerical 
complexity of the problem. Moreover, we may want to 
extend our basis set beyond the nine bands 4s, 4p, and 
3d to allow for a further relaxation of the one-particle 
wave functions. Again, this is feasible but numerically 
much more costly. Lastly, one may want to include some 
atomic p-p and p-d interactions. All these extensions are 
left for future studies. 

C. Comparison to experiments 

Typical results of our calculations for Ni are summa- 
rized in table [I]. They are obtained from the multi- 
band Gutzwiller method utilizing the nine-orbital tight- 
binding model based on DFT energy band calculations 
for non-magnetic Ni, and employing values of A = 10 eV 
and C = 0.1 eV, with C/B = 4.5. 

The width of the d bands is predominantly determined 
by A (essentially the Hubbard-C/), via the values of the 
hopping reduction factors qb, a - The exchange splittings 
and, consequently, the magnetic moment are strongly in- 
fluenced by C and only moderately by A. The Racah 
parameter C causes the Hund's-rule splitting of the d s 
multiplets; in the hole picture, d s is the only many- 
particle configuration which is significantly occupied (by 
1.90 electrons), while 5.94 electrons are in d 9 , 0.89 elec- 
trons are in d 10 , and 1.18 electrons have s or p character. 

Generally, the Gutzwiller results agree much better 
with experiment than the SDFT results. This is the case 
for, (i), the d bandwidth W* , (ii), the value for N*(E F ), 

(iii) , the positions of individual quasi-particle energies, 

(iv) , the values of the exchange splittings, (v), their ti g - 
e g anisotropy, and, (vi), the ti g je g ratio of the d part 
of the magnetic moment. As a consequence of the small 
d(e g ) exchange splitting, the X21 state lies below E-p and, 
thus, the Fermi surface exhibits only one hole ellipsoid 
around X, in nice agreement with experiment. 



G 





Exp 


Ref. 


this work 


SDFT 


/•^spin-only / /-^B 


0.55 | 


U 




0.53 


0.59 


r(t 2g /e s ) 


4.3 


1 


34 




5.1 


2.8 


W* cc = -(Xi)/eV 


-3.3 ± 


0.2 [ 


2£ 




-3.14 


-4.86 


W*/eV 


3.5 ± 


0.3 [ 


2£ 




3.29 


5.20 


N*{E F )/(eV atom) 


3.0 


1 


3C 




2.6 


1.9 


(X 3 )/eV 


-2.8 ± 


0.2 [ 


K 




-2.94 


-4.06 


X 2T /eV 


9/1-1- 


n r»9 I 


42 


.55 


-0.22 


-0.56 


X 2l /eV 


-0.06 ± 


0.02 


1 


,3E 


-0.04 


+0.09 


X 5T /eV 


-0.10 ± 


0.02 


i 


,44 


; -0.14 


-0.22 


Xsi/eV 


0.23 ± 


0.04 




,44 


; 0.25 


0.34 


<r 12 >/ev 


-0.4 ± 


0.1 






-0.64 


-1.02 


<r 25 '>/eV 


-1.1 ± 


0.2 






-1.45 


-2.14 


<L 3 )/eV 


-1.3 ± 


0.1 






-1.49 


-2.23 


A((l/2)ri)/eV 


0.21 ± 


0.02 






0.30 


0.60 

















TABLE I. Comparison of experimental and theoretical re- 
sults. Parameters for the present work are: A = 10 eV, 
C = 0.1 eV, and C/B = 4.5. The SDFT results are 
quoted from Ref. [^J; the result for r(t2 g /e g ) is taken from 
Ref. j}9|. /ispin-oniy: spin contribution to the magnetic 
moment; r(i2 S /e s ): ratio between t 2g and e g contribution 
to the magnetic moment pi S pin-oni y ; W* = X51 — (Xi): 
quasi-particle bandwidth; W^ cc : occupied quasi-particle 
bandwidth; N*(Ef): quasi-particle density of states at _Ef; 
(. ..)/eV: spin-averaged energies of the quasi-particle states 
at various high-symmetry points of the Brillouin zone; X 2i 5 CT : 
energies of the spin-a quasi-particle bands with respective 
symmetries X2/-X5 at the X point; A((1/2)FL): typical ex- 
change splitting of the highest-lying A3 quasi-particle states 
halfway between V and L. 

The large anisotropy of the exchange splittings is a 
result of our ground-state energy optimization, which 
allows A f2g and A eg to be independent variational pa- 
rameters. We find A t2 „ w 3A e „ w 800 meV. Note 



r 2g 

that these values enter |<I>Q arc ) and are renormalized by 
factors qi,^ w 0.7, (ftj ~ 0.6, when |$Q Pt ) is reached. 
This also implies that the width of the majority spin 
bands is about 10% bigger than that of the (higher ly- 
ing) minority spin bands. It causes a further reduction 
of the exchange splittings of states near Ep, especially 
for those with strong t 2g character. Note that this band 
dispersion effect causes larger exchange splittings near 
the bottom of the d bands, e.g., 0.45 eV splitting of X\ 
and 0.74 eV splitting of X 3 . There, however, the quasi- 
particle linewidths have increased to 1.25 eV and 1.4 eV, 
respectively 29 1, so that an exchange splitting near the 
bottom of the d bands could, so far, not be observed 
experimentally. When we investigate the quasi-particle 
energies in the middle of the d bands, such as Ti 2 , r 2 5', 
or L3, we observe that our theoretical values lie consis- 
tently below the experimental ones. If we assume that 
these states can be measured more reliably than those 
near the bottom of the d band, we may conclude that 
the value of A could still be somewhat bigger than 10 eV. 
Studies with larger A values are presently carried out. 

We think that the large anisotropy originates from pe- 



est neighbors, and 



culiarities special to Ni with its almost completely filled 
d bands and its fee lattice structure. Near the top of 
the d bands, the t 2g states dominate because they ex- 
hibit the biggest hopping integrals to nearest neighbors, 
The e„ states have t^L ~ —0.3 eV to near- 
dda ~ 0.1 eV to next-nearest neigh- 
bors; the latter are small because of the large lattice dis- 
tance to second neighbors. The e g states also mix with 
the nearest-neighbor t 2g states with t^ n -type coupling. 
Therefore, the system can gain more band energy by 
avoiding occupation of anti-bonding t 2g states in the mi- 
nority spin bands via large values of At 2 , at the expense 
of allowing occupation of less anti-bonding e g states via 
small A e values. Note that this scenario is special to 
nickel and should not apply to materials with a bec lat- 
tice structure which have almost equal nearest and next- 
nearest neighbor separations. Since the bands in nickel 
are almost completely filled, the suppression of charge 
fluctuations actually reduces the number of atomic con- 
figurations where the Hund's-rule coupling is active. It 
is also in this respect that nickel does not quite reflect 
the generic situation of other transition metals with less 



completely filled d bands; see Sect. IV. 

The agreement between theory and experiments on 
nickel shows that our multi-band Gutzwiller theory is a 
useful tool for the theoretical description of the ground- 
state properties of transition metals and their com- 
pounds. In the next section we describe a way to calcu- 
late consistently low-energy excitations in itinerant fer- 
romagnets within our variational approach. 



VI. SPINWAVES IN ITINERANT 
FERROMAGNETS 

As in magnetic insulators, the elementary excita- 
tions in itinerant ferromagnets are spinwaves. Since 
ground-state properties are well described by Gutzwiller- 
correlated wave functions it is highly desirable to derive 
spinwave dispersions from the same variational approach. 

In fact, the variational principle can also be used to 
calculate excited states Eqj. If |$) is the ferromagnetic, 
exact ground state with energy Eq, the trial states 



(19) 



are necessarily orthogonal to |<E>), and provide an exact 
upper bound to the first excited state with momentum q 
and energy e(q) 



< E s (q) 



_ (<f( q )\H\*(q)) 



En 



(20) 



Here, S q = (S+)+ = J2i,b ex P( -i< ?0C&,i^,&,T flips a spin 



from up to down in the system whereby it changes the 
total momentum of the system by q. In this way, the 
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famous Bijl-Feynman formula for the phonon-roton dis- 
persion in superfiuid Helium was derived [ [46[ . In the case 
of ferromagnetism the excitation energies E s (q) can be 
identified with the spinwave dispersion if a well-defined 
spinwave exists at all p7[ . Experimentally this criterion 
is fulfilled for small momenta q and energies E s (q). 

Unfortunately, we do not know the exact ground state 
or its energy in general. However, we may hope that the 
Gutzwiller wave function |^g) is a good approximation 
to the true ground state. Then, the states 



|* G («)) =S-\* Q ) 
will provide a reliable estimate for E s (q), 



(21) 



E s (q) a El 



= (^ G \S+HS-\^ G ) 
(* G |5+4-|* G ) 



(gGjg|gg) 

(*g|*g> 



(22) 



Naturally, E™ (q) does not obey any strict upper-bound 
principles. 

The Hamiltonian (Q) commutes with the operator for 
the total spin. Thus, for q = 0, S~_ Q generates a state 
from the same spin multiplct. Therefore, the spinwave 
at zero momentum is the Goldstone mode of the ferro- 
magnetic system, E s (q = 0) = 0. Within the variational 
approach E™'(q — 0) = is guaranteed if we choose 
I^g) as an eigenstate of the total spin operator. This 
is most easily accomplished by a minor restriction of the 
variational space for the atomic occupancies J27j with 
negligible influence on the variational ground-state prop- 
erties. 

The actual calculation of the variational spinwave dis- 
persion is rather involved even in the limit of large di- 
mensions. We note, however, that explicit formulae are 
available |27|] which can directly be applied once the vari- 
ational parameters have been determined from the min- 
imization of the variational ground-state energy. As a 
first application, we present results for the two-band toy- 



model of Sect. IV. 

In Fig. we show the variational spinwave dispersion 
in x direction, E^ aT ((q x ,0,0)), for the model parame- 
ters n a — 0.29, J = 0.2/7, and the four different val- 
ues U/eV — 7.8, 10, 12, 13.6 which correspond to a mag- 
netization per band of m = 0.12,0.20,0.26,0.28. This 
quantity is defined as < m — (nt,^ — n&.|)/2 < n/4. 
Note that our last case corresponds to an almost com- 
plete ferromagnetic polarization. The data fit very well 
the formula 



E™((q x ,0, 0)) = Dql(l + Pql) + 0(q x ) , 



(23) 



in qualitative agreement with experiments on nickel E7J . 
The corresponding values D = 1.4 eVA 2 and D = 
1.2 eV A 2 for m = 0.26 and m = 0.28, respectively, 
are of the right order of magnitude for nickel where 
D = 0.43 eVA 2 . As lattice constant of our simple-cubic 
lattice we chose a = 2.5A. 




q/A" 

FIG. 3. Variational spinwave dispersion in x direct ion, 



£r r ((<2z,0,0)), for the two-band model denned in Sect. |V 
n/4 = 0.29, J = 0.2(7, and the values U/eV = 7.8, 10, 12, 13.6 
correspond to m — 0.12, 0.20, 0.26, 0.28. The lattice constant 
is a = 2.5A. Inset: £ s var ((g*, 0, 0)) and £ s var ((g*, q x , 0)) for 
m = 0.2 and m = 0.28, respectively. The spinwave dispersion 
is almost isotropic. 

As shown in the inset of Fig. the dispersion relation 
is almost isotropic for q x values up to half the Brillouin 
zone boundary |27j ], in particular for large magnetiza- 
tions. This is in contrast to the strong dependence of the 
electron-transfer amplitudes tjj on the lattice direction. 
This implies for strong ferromagnets that the collective 
motion of the local moments is similar to that of localized 
spins in an insulator ply . Such ferromagnetic insulators 
are conveniently described by the Hciscnbcrg model with 
exchange-interaction between neighboring sites on 
a cubic lattice, 



(24) 



(i,3) 



For such a model one finds D = 2SJa 2 . The length 
of the effective local spins can be calculated from |Wq) 
as S{S + 1) w 0.95 (S = 0.6) for m > 0.20 @. 
Therefore, J s» D/(1.2a 2 ), which gives the typical value 
J = 0.17 eV. For an estimate of the Curie tempera- 
ture Tc we use the result from quantum Monte-Carlo 
calculations E9I 



T c = 1. 44 JS 2 



(25) 



for spins 5 on a simple-cubic lattice. In this way we find 
T c ~ 0.5 J = 0.09 eV = 1 • 10 3 K. This is the same or- 
der of magnitude as the condensation energy for these 
values of the interaction, i? C ond = 5 • 10 2 K, see Sect. IV. 
Given the arbitrariness in the relation between -E C ond and 
Tc, and the application of the Heisenberg model to our 
itinerant-electron system, we may certainly allow for dif- 
ference of a factor two in these quantities. Nevertheless, 
the results of this section clearly show that, (i), -E C ond 
gives the right order of magnitude for Tc, and that, (ii), 
the spinwave dispersion of strong itinerant ferromagnets 
resemble the physics of localized spins. 
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VII. CONCLUSIONS 

Which scenario for itinerant ferromagnetism in transi- 
tion metals is the correct one? 

For a long time this questions could not be answered 
conclusively. Band theory along the lines of Slater and 
Stoner could be worked out in much detail whereas a 
correlated-electron description of narrow-band systems 
was lacking until recently. In fact, van Vleck stated in the 
summary to his conference contribution in 1953: "The 
gist of this paper is that it would be highly desirable if 
good methods of computing with (c) (minimum polarity 
model) could be developed." In 1963, Gutzwiller laid the 
foundations for a concise treatment of the problem, and 
today we are able to draw definite conclusions based on 
Gutzwiller-correlated multi-band wave functions. 

Our results for nickel indicate that the van-Vleck- 
Gutzwiller scenario is valid. Band theory does not ap- 
preciate the strong electronic correlations present in the 
material which lead to the observed renormalization of 
the effective mass, exchange splittings, bandwidths, and 
Fermi surface topology. Moreover, as also covered by 
our approach, charge fluctuations are small, and large 
local moments are present both in the paramagnetic and 
the ferromagnetic phases. Roughly we may say that the 
electrons' motion through the crystal leads to a ferromag- 
netic coupling of pre-formed moments which eventually 
order at low enough temperatures. In this way, strong 
itinerant ferromagnets resemble ferromagnetic insulators 
as far as their low-energy properties are concerned: spin- 
waves exist which destroy the magnetic long-range order 
at the Curie temperature. 

It is thus seen that today computational difficul- 
ties no longer obscure "the recognition in principle of 
the situation which confirms closest to physical real- 
ity" (van Vleck, 1953 ||): transition metals and their 
compounds are strongly correlated electron systems. 
More than thirty-five years after their introduction by 
Gutzwiller 0, our studies clearly show that Gutzwiller- 
correlated multi-band wave functions successfully de- 
scribe the low-energy physics of these materials. 
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